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ABSTRACT 

Late-time spectra of stripped-envelope CC-SNe are dominated by strong [O i] 
AA6300,6363 emission, caused by thermal electron excitation of forbidden [O i] tran- 
sitions. The permitted O i 7774 A line is also often observed. This line cannot result 
from thermal electron excitation of the oxygen ground state. In this work tests are 
performed to verify whether the line can be powered by oxygen recombination alone, 
using the examples of two of the best studied type Ic SNe, 1998bw and 2002ap. 

Temperature-dependent effective recombination coefficients for neutral oxygen are 
calculated using available atomic data. Missing atomic data are computed in a temper- 
ature range typical for SN nebulae. Core ejecta models for SNe 1998bw and 2002ap are 
obtained from modelling their nebular emission spectra so that oxygen recombination 
line formation is computed consistently with oxygen forbidden line emission. 

While SN 2002ap can be explained well by a one dimensional shell model, this 
seems not to be possible for SN 1998bw, for which a two dimensional model is found. At 
very late epochs the formation of the O I 7774 A line can be explained by recombination 
radiation for both SNe, but at earlier epochs strong absorption is present which may 
determine the strength of this line even at ~ 200 days. 
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I INTRODUCTION 

Massive stars (> 8 Mq) collapse when the nuclear fuel in 
their central regions is consumed, producing a core-collapse 
supernova (CC-SN) and forming a black hole or a neutron 
st ar. CC- SNe with a H-rich spectrum are classified as Type 

II (|Filippen ko 1997,). If the star was stripped of at least most 
of its H envelope prior to the explosion, the SNe are classified 
according to the degree of stripping as Type lib (strong He 
lines, and weak but clear H), Type lb (strong He lines but 
no H), and Type Ic (no He or H lines). 

Some CC-SNe, called broad- lined SNe (BL-SNe), ex- 
hibit very broad absorption lines at early times, caused 
by the presence of sufficiently massive ejecta expanding 
at high velocities. These SNe are sometimes called hyper- 
novae, and can be a ssociated with long-durat ion gamma-ray 
bursts (GRBs) (see lWooslev fc Bloom|[200^ . and references 
therein). 

In GRB scenarios a relativistic outflow launched by a 
central engine deposits some fraction of its energy into the 
SN ejecta. Since the energy is probably deposited preferen- 
tially along the polar axis, the SN might be strongly asym- 



metric (e.g.. iMaeda et ai]|2002l ). The neares t, best-studied 
GRB -SNe are SN 1998bw / GRB 980425 (ICalama et all 
|l99i), SN 2003dh / GR B 030329 jMathesonl l2004h . SN 

' 112004'), and SN 

I2OO6I . 



20031w / GRB 031203 (' Malesani et al 
2006aj / GRB/XRF 060218 jPian et al 



yet fully established whether the GRBs (or X-ray flashes) 
accompanying nearby CC-SNe share the same properties 
of high-redshift GRBs. CC-SNe may also be characterised 
by asphericities, although a jet does not necessarily form 
(iBlondin et al.l l2003l: iKotake et all |2004 fMoiseenko et all 



I2OO6I : iBurrows et alj|2007l : iTakiwaki et al.ll2009l ) 



Asphericities in the inner and outer ejecta are evident in 
at least some CC-SNe. Two clear indicators a re velocity dif- 
feren ces of Fe and lighter-element lines (e.g.. iMazzali et al.l 
l200lh . and polarisation measurements (e.g.. IlKfiichl 1991 ). 
Indirect indications also emerge from a comparison of the 
inner and outer ejecta velocities (jMaurer et al.1120091 ). 

Independent of their type, SNe become increasingly 
transparent to optical light with time, as the ejecta thin out. 
At late times (> 200 days after the explosion), the inner- 
most layers of the SN can be observed. This epoch is called 



© 0000 RAS 



2 Maurer et al. 



the nebular phase, because the spectrum turns from being 
dominated by absorption to an emission spectrum, mostly 
showing forbidden lines. In this phase the radiated energy of 
a SN is provided by the decay of radioactive ^^Co (which is 
produced by the earlier decay of ^®Ni) into ^®Fe. Decaying 
^®Co emits 7-rays and positrons which are absorbed by the 
SN ejecta. As the deposition rate of 7-rays and positrons de- 
pends on the density and *^Ni distribution, the inner parts of 
the SN dominate the nebular spectra. Therefore the nebular 
phase is especially suitable for studying the core of SNe. 



Several authors have modelled nebular-phase spec- 
tra of SNe to derive q uanti ti es such as the ^^N i 



mass (e.g.. 


Mazzali et al.l 12004 


: IStritzineer et all I2OO6I: 


Sauer et al. 


20061: iMaeda et al 


2007 


d), ejecta veloci- 


ties fe.e.. iMazzah et al.l l2007bl: 


Taubenbereer et al.l 


200g|: 


Maurer et al. 


I2OO9I'). asDhericities fe.e.. Mazzali et al.l 


2005I: 


Maeda et al.ll2006l: iMazzali et al. 


2007al. 


200a). and elemen- 


tal abundances fe.e.. iMaeda et a" 


.112007a 


Ibt). 





To infer SN core ejecta velocities and geometry usually 
the [O l] AA 6300, 6364 doublet is investigated, since this 
is by far the strongest nebular emission line (SNe of type 
Ic), especially at very late epochs. These lines are formed 
by thermal electron excitation of the 2p(^D) state. Higher 
quantum states (n 3) cannot be excited by thermal elec- 
trons sufficiently because of the large ratio of excitation en- 
ergy (~ 9 eV) and electron temperature (~ 0.4 eV). Possi- 
ble excitation mechanisms are absorption, non-thermal elec- 
trons and recombination. 

Previous studies of central oxygen have focused almost 
exclusively on the [O i] AA6300, 6364 doublet. It is useful to 
investigate whether extending the analysis to other oxygen 
lines gives a consistent picture. Therefore it is important to 
understand the formation of these lines in detail. 

In Section 2 temperature-dependent effective recombi- 
nation rates are obtained for the n = 3 levels of neutral 
oxygen and some quantities relevant for oxygen recombi- 
nation and emission-line formation in CC-SN nebulae are 
estimated. A one dimensional shell model for SN 2002ap is 
described in Section 3 and a two-dimensional model for SN 
1998bw in Section 4 emphasising the role of the O I 7774 A 
line. Results are discussed in Section 5. 



2 OXYGEN LINES IN THE NEBULAR PHASE 
OF CC-SNE 

In this section temperature-dependent effective recombina- 
tion rates are obtained for the 3s, 3p and 3d triplet and 
quintet states of neutral oxygen. These rates are used in our 
nebular code to calculate oxygen recombination line emis- 
sion. These levels are responsible for O I 7774, 8446, 9264 
and 11287 A line emission (see Figure[T]for illustration). We 
further investigate the relevant processes for oxygen recom- 
bination and absorption line formation in the nebular phase 
of CC-SNe. 

Recombination into the 2p singlet states is not impor- 
tant in the CC-SN nebular phase. About 80 - 90% of the 
radioactive energy is deposited in thermal electrons, while 
only 10 — 20% goes into ionisation (see below). Since the 2p 
singlet states can be efficiently excited by thermal electron 
collisions, the population of these levels is predominantly 
determined by this process. 




Figure 1. Neutral oxygen levels that are potentially interesting 
for the nebular phase of CC-SNe. The singlet states (which are 
responsible for the nebular oxygen forbidden-line emission) and 
states with n ^ 4 are not shown. The 2p, 3p and 3d states can 
be subdivided into several j sub-states, which are not shown to 
keep the illustration clear. The wavelengths given represent the 
mean wavelengths of emission. Lines emitted by the different j 
sub-states differ slightly from the values given. For example, the 
line we refer to as O I 7774 A in this paper consists of three 
lines at wavelengths 7771.94, 7774.17 and 7775.39 A. While all 
transitions indicated by solid lines are permitted with radiative 
rates of order 10^~* , the 3s(^S) to 2p('^P) transition (spin- 
forbidden; dashed line) has a radiative rate of order 10'' only. 



2.1 Effective recombination rates for neutral 
oxygen 

Di-electronic recombination is much weaker than radiative 
recom bination [roughly a factor of 5 (Nussbaum er fc StorevI 
11983) at 10000 K and decreasing rapidly with tempera- 
ture] and is not included when calculating the recombina- 
tion fractions. At high electron densities coUisional transi- 
tions between excited states can become important. How- 
ever, for the electron densities expected in the nebular 
phase (tIc < 10^'^ cm~^), this effect is weak. At high elec- 
tron densities (no ~ 10^" cm"'') and low temperatures 
(r ^ 3000K) coUisional recom bination can become impor- 
tant (e.g. see lStorev fc Humm er 1995). However, during the 
nebular phase only the combinations high density, high tem- 
perature (early) and low density, low temperature (late) are 
important and the contribution of coUisional recombination 
is weak (e.g. ~ 20% at ~ 10^ cm~^ and T ~ 5000K). 

Therefore we only consider radiative recombination. We 
are interested in which fraction of recombining electrons 
reaches a certain atomic level, whether by direct recombina- 
tion or cascad i ng fro m higher state s . It h ad been noted by 
lJulienne et~aLl l| 19741 ') : [Chung et al.l ([l99j) that high quan- 
tum levels (n ~ 20) might be important to calculate the 
recombination cascade. In the literature recombination co- 
efficients and radiative rates are available for neutral oxygen 
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large radii and using a Thomas-Fermi potential 
r r 



(2) 
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Figure 2. Temperature-dependent normalised effective recom- 
bination fractions. The quintet states (3d, 3p, 3s from bottom 
to top) are shown by solid lines. Triplet states (3p, 3s, 3d from 
bottom to top; note that 3p and 3s do overlap) for the case of 
optically thin ground state transitions are shown by dashed lines. 
Triplet states (3d, 3p, 3s from bottom to top) for the case of op- 
tically thick ground state transitions are shown by dotted lines. 



up to maximum quantum numbers n ~ 10 (see e.g. TIPTOP- 
baseB NIStQ). 

For azimuthal quantum numbers 1^3 one can use the 
hydrogenic approximation, howe ver for the s, p a nd d or- 
bitals this should be avoided (e.g. lChung et alll99ll ). There- 
fore, we obtained radiative recombination coefficients and 
radiative rates for all s, p and d levels of oxygen up to n = 20 
using a quantum d e fect method (QDM) |Bates fc DamgaardI 
1 19491 : ISeatonll 19581 : iBurgess fc Seatonlll960h . 

The potential of the atomic core and the inner electrons 
approaches the asymptotic form 1 /r quickly. The wave func- 
tions of the excited states can be described by a hydrogen- 
like solution, with the difference that the quantum energy 
levels are shifted owing to the unknown structure of the 
core potential. This shift of energy levels is described by the 
quantum defect (defined as the difference between the real 
and 'effective' principal quantum number) and can be mea- 
sured in experiments or is determined by more sophisticated 
calculations. 

The wave functions calculated by the QDM are de- 
scribed with sufficient accuracy at radii larger than a few 
(atomic units) and therefore the QDM is reliable as long as 
the main contribution of the transition integrals comes from 
sufficiently large radii and there is no strong cancellation in 
the integrand. Especially for ground state transitions, where 
contributions from small radii are important, the QDM is 
therefore not reliable. 

To improve the calculation , we numeri cally integrate 
the electron wave equation (e.g. ISeatoniri958l ) 



^(L^^V{r) + E\P{E,l-r)=0 (1) 



inwards, starting with the solution obtained by the QDM at 

^ http://cdsweb.u-strasbg.fr/topbase/ 
^ http://www.nist.gov/indcx.html 



where Z is the charge of the nucleus and N the number 
of core electrons. This integration is not part of the QDM. 
The Thomas-Fermi potential gives the correct asymptotic 
beha viour of the atom ic potential for r — >■ and r — >■ oo 
( Bur ke fc RobbI Il975l ). The result of this integration is a 
behaviour of the wave-function for r — >■ different than pre- 
dicted by the QDM. For most transitions the contribution 
from small radii is approximately zero and therefore the 
calculated rates are not inffuenced by this correction. For 
the transitions which are inffuenced by the inner parts of 
the wave function we achieve better agreement with oxy- 
gen atomic data provided in the literature at low quantum 
numbers (n < 10). 

The difference for ground state radiative transitions rel- 
ative to NIST recommended data is always less than 50% 
and the agreement increases considerably with the quan- 
tum number of the excited states. For ground state transi- 
tions from levels n > 5 the disagreement becomes less than 
10%. For transitions between excited states the maximum 
difference to NIST recommended data is 40% with most 
transitions agreeing at the 10% level or better. Whenever 
the disagreement was worse than 10% we replaced our data 
with the values from the literature. We comp ared our direct 
and effective recombinatio n coefficients to ljulienne et al.l 
(|l974h : IChung et al.ll|l99ll) and found good agreement (bet- 
ter than 10%), as we did for the effective recombination 
rates at 1160 K. Our total re combination rate agrees with 
the value given bv lChung etal . (1991) to 1 % at 10000 K. 
However compared to lAldrovandi fc Peguignotl (|l973l ) our 
tot al recombination rat e is ~ 25% too low. It was noted 
bv I Julienne et al] (|l974h that neglecting the n > 20 levels 
might cause an underestimate of ~ 20% of the total recom- 
bination rate (at temperatures around 1000 K; the effect 
becomes smaller at highe r temp eratures). More importantly 
lAldrovandi fc Peguignotl l| 19731 ') used a hydrogenic approxi- 
mation for all atomic state s, which leads to an overestimate 
of the recombination rate (IChung et al.|[l99ll). A compari- 
son to a recent calculation from Badnel 1 (|2006t ) shows that 
our total radiative recombination rates agree to their results 
within 15% . To handle that deviation, we normalise our re- 
sult. 

We calculate the ratio of the effective recombination 
rates of the n = 3 levels to the total recombination rate into 
excited states and find that this ratio varies only weakly with 
the number of levels included (it changes by ~ 1% (3s), 1% 
(3p) and 8% (3d) when using a maximum quantum number 
of n = 20 instead of n = 10) and is not too sensitive to 
uncertainties in the atomic data (see below), while the total 
recombination rate is (it increases by 15 - 30% depending 
on temperature when using a maximum quantum number n 
= 20 instead of n = 10). 

Since fitting formulae for the total and ground state 
recombination rates are available in the literature, we can 
obtain the total recombination rate into excited states and 
normalise our effective recombination rates to the total ra- 
diative recombination rate. We call this the normalised ef- 
fective recombination fraction. 

The recombination rates used for this normalisation are 
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Figure 3. The ratio fio[l]/'^Tot for different clumping factors 
= 1, 10, 100 from bottom to top) calculated for the test model 
described in Appendix A. Since clumping increases recombination 
the fraction of neutral oxygen increases with C,. 



the total recombination rate of iBadnell (|2006l l 



-Ro[il — ^ 



(3) 



-11 



B = 



A = 6.622 ■ 10" 
4.136, Ti = 4.214 ■ 10^ and 



with B' = B + Cexp ( 
0.6109, C = 0.4093, To 
T2 = 8 .770- 10'* and the (direct) ground state recombination 
rate of lPeauigno^ (|l990l ') 



Rgs = (1.174 + 0.2463a; + 0.2144a;^ - 0.0621a;^ 



w in-13rrn-0.5 3 -1 

X 10 Ta cm s 



(4) 



with X = log(l+r4)/ log(2) and T4, = T/10000 K. These nor- 
mahsed effective recombination fractions are shown in Fig- 
ure [2] and are hsted in Appendix B. To obtain absolute val- 
ues one has to multiply by the total radiative recombination 
rate. Since the total recombination rate contains a contribu- 
tion from direct ground state recombination and since lower 
levels contain cascading contributions from higher states, 
these fractions do not add up to one. 

Since there are uncertainties in our atomic data we want 
to quantify their influence. We randomly vary all the radia- 
tive rates obtained by the QDM by a value between ± 20% 
and the recombination rates obtained by the QDM by a 
value between ± 50% several times. This should be larger 
than the typical error, at least for excited states. The result- 
ing deviations from the standard values are typically a few 
percent, often <C 5%. 

2.2 Recombination line formation 

With the effective recombination rates shown in Figure [2] it 
is possible to calculate the strength of all [O l] n = 3 lines 
due to recombining electrons. 

The population of excited oxygen levels by recombina- 
tion i s imple mented in the one dimensional IIMazzali et al 



2001 , l2007al ) and three dimensional version (Maurer et al. 



2009|) of our nebular code and will be used in Sections [3] and 



4]to calculate synthetic recombination lines for SNe 1998bw 
and 2002ap. 



However, to obtain insight into the formation of these 
lines it is useful to derive some estimates for the recombi- 
nation line formation. Below, we estimate the strength of 
the O I 7774 A line. This is the strongest non-ground state 
(NGS) recombination line thanks to its high effective re- 
combination rate and transition energy [3p('''P) to 3s(^S)]. 
Ground state recombination line observations are usually 
not available since most spectra reach a minimum wave- 
length of 3000 A - 4000 A only. 

The ionisation rat e can be c alculated using the concept 
of 'work per ion ' (e.g. lAxelrod|[l98Q ). which describes the 
energy lost to thermal electrons in order to produce one 
ionisation. The ionisation rate per atom and ion is then given 
by 



Yo 



o[i] 



(5) 



where Ldgp is the total deposited luminosity, A^xot is the to- 
tal number of atoms and Wo[/] is the energy lost to thermal 
electrons per ionisation and can be estimated by compar- 
ing the ionisation cross-section with atomic and plasma loss 
functions (more details are given in Section [23} . 

During the first several hundred days of a SN, ionisation 
balance is a valid assumption (Axolrod 1980!), therefore 

Yo[i](l + ■^)'^o[Il = io[i]'^o[i] = -Ro[I]^^cC'^o[IIl (6) 

where TZ is the recycling fraction (recombination radiation 
causes further ionisation and TZ is expected to have a value 
between 0.3 and 0.5; lAxe lrod 1980) and ^ is the clumping 
factor (explained in mor e detail below), defi ned as the in- 
verse of the filling factor (|Li fc McCravlll993 ). 

Since oxygen is the dominant element in CC-SNe cores 
[~ 55 % in SN 2002ap and - 75 % in SN 1998bw, as ob- 
tained from our modelling in Sections [3] and [4] see also 
iMazzali et"aLl (|200ll . l2007al )] and the degree of ionisation 
of oxygen is rather high (roughly 10 — 50%; Figure[3]) it is a 
reasonable simplification to neglect all elements besides the 
iron-group elements and oxygen when calculating the ioni- 
sation balance and to assume that all electrons are provided 
by oxygen and by iron, which is assumed to be 100% singly 
ionised. Therefore 



yO[I]'^0[I] ~ -Ro[I]('^0[II] + 'lFc)C"-0[IIl 



(7) 



which gives 
nom , 



riTot 



npc Ro[i]nFcC + Yq 
riTot 



o[i] 



2i?o[I]WTotC 



1 + 



^ ^ 4y(^[i]-Ro[i]('^T, 

(i?o[i]»^FcC + i^o[r 



- npc)C^ 



n-Tot 



-Ro[i]nFcC + yo[i] 
2i?o[I]»^TotC 



y"o[i]("Tot — npo 



-RoiiiJ^TotC 
riTot 

y(nTot — npey 
(i?0[I]'^FeC + yo[I])riTot 



"^o[il('^Tot — npc 
-Ro[i]w|otC 



> 1 



yo[i](wTot - npc) 

-Ro[I]"-TotC 



< 1 



(8) 

In Figure |3]this ratio is computed for a test model described 
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in Appendix A using clumping factors of 1, 10 and 100. 
The total luminosity of any recombination line X of oxygen 
relative to the total deposited luminosity is then given by 

ix 



^Dcp 



= io[I]'JO[I]V'Sx/xI/Dcp 



= (l + 7^) 



Ex fx fiT'Oii] 



Worn V TiTot 



(9) 



where Ex is the energy and fx is the normahsed effective 
recombination fraction of the line. Usually one can assume 
4:Yni^Rn mnTotC > (fin [T]nFeC + ^"^[7])^ ; see Appendix A. 
Loep fe.g. lAxelrodlll98ol ) and all other variables (besides the 
temperature) can be calculated for any ejecta model and 
Equation [9] can be evaluated directly. The temperature can 
be determined by balancing heating and cooling taking into 
account several hundred emission lines (therefore an exact 
estimate of the temperature is difficult) and varies between 
8000 K and 2000 K at the epochs of interest (100 to 600 
days). A sufficiently exact estimate can be obtained assum- 
ing a constant temperature of ~ 4000 K (see Appendix A). 

The luminosity of any oxygen recombination line is di- 
rectly proportional to the normalised effective recombina- 
tion fraction fx and the energy of the transition Ex [eV]. 
The recombination line luminosity is influenced only weakly 
by the clumping factor (which changes the ionisation frac- 
tion). The oxygen NGS recombination lines will be weak in 
general [taking no[i]/nTot ~ 0.5, /7774 ~ 0.4, £7774 ~ 1.6 eV 
and Wo [I] ~ 75 eV the strongest NGS oxygen recombination 
line (7774 A) carries less than 1.0% of the total luminosity] 
and we will neglect all other NGS oxygen recombination 
lines in the rest of this paper. They are calculated by the 
nebular code, but in the nebular spectra of SNe 1998bw and 
2002ap on oxygen recombination lines other than 7774 A 
can be identified. 



2.3 Excitation of the O I 7774 A line 

We will find in Sections[3]and|4]that the luminosity provided 
by recombination is not sufficient to explain the observations 
of the 7774 A line earlier than ~ 250 days. Therefore, in this 
section we investigate other excitation mechanisms which 
might operate simultaneously with recombination. 

Ground-state excitation by thermal electrons can 
clearly be ruled out. All transitions from n 3 to the ground 
state have energies larger than 9 eV, which makes thermal 
excitation by electrons of temperatures ~ 0.5 eV ineffective. 

Non-thermal excitation by the same electrons causing 
ionisation (which are produced by Compton-scattering of 7- 
rays emitted by ^''Co decay) is also too weak. The ratio of the 
cross-sections for electron impact excitation and ionisation 
(at the high electron energies which arc of interest here) can 
be approximated by (e.g. Rozsn vai ot al...l980. ) 

^ni^ ^ 1.66if^H^ (10) 

where E^i^nn' is the transition energy, Eni the ionisation 
energy, A^ni is the number of electrons in the shell n, 1 and 
fni-m'i' the oscillator strength of the transition. We compare 
hydrogen and helium non-thermal excitation rates obtained 
by this approximation to values obtained from more sophis- 
ticated calculations (Hachinger in prep.), solving an energy- 



balance equation der i ved from t he Spencer-Fano equation 
(|Xu fc McCravlll99ll : lLucvlll99lh and find that the agree- 
ment is always better than 30%, which suggests that this 
approximation is also of acceptable accuracy for oxygen. For 
allowed ground state transitions of neutral oxygen the os- 
cillator strengths are of the order of 10~^ decrea sing with 
mcreasmg quantum number (|Bell fc Hibbertlll990l ). The to- 
tal excitation rate into triplet states is ~ 5% of the ioni- 
sation rate only. The excitation rates of quinted states by 
non-thermal (high-energetic) electrons is even lower, since 
the excitation cross-sections of spin-forbidden lines strongly 
decrease wi th electron energy, as c ompared to allowed tran- 
sitions (e.g. iRalchenko et al.ll2008l ). 

Therefore, non-thermal excitation is not important 
compared to recombination for populating excited levels. In 
addition, if the O I 7774 A line was excited by non-thermal 
electrons, one would expect the ratio of the line to the to- 
tal luminosity to be approximately constant at all epochs, 
which is clearly not observed (this ratio in fact decreases 
considerably at late epochs). 

Although it is generally assumed that the SN is optically 
thin during the nebular phase, some lines can still be opti- 
cally thick (especially ground state transitions). However, as 
we show here, the optical depth of the 7774 A transition is 
still high, even at epochs of ~ 200 days because the 3s(^S) 
level is populated by recombining electrons. 

We estimate the time-dependent optical depth of the 
7774 A line resulting from recombination taking into account 
the effective recombination rates, resonance scattering of the 
3s(^S) to 2p(''P) transition as well as the recycling of ground 
state recombination radiation. 

Based on Equations [7] and [8l the number density of the 
3s(^S) state relative to the total density is approximately 
given by 

""■SsC^S) /3s(5S)-Ro[I](wFc + Mo[ii])Cwo[IIl 



nTot (A' + C)nTot 

A V nxot 



v4 > A' > C 



(11) 



(l + 7^) 



ylV"Wo[Il V JlTot 



defining A' as the reduced (due to self-absorption) radiative 
rate 



,, , 1 — exp(— r) A 
A = A — ~ — , r > 1 



(12) 



from the 3s(^S) to the ground state, where r is the Sobolev 
optical depth of this ground state transition 

-^1355 *S3s{5 3)^1355^0 [I] C 



87r52p(3p) 



To(no[i] , no[i] » n-3s(5s) 
(13) 

with the ground state transition wavelength A1355 = 1355 A, 
t the epoch in seconds, V the volume in cm^, A1355 the ra- 
diative rates from the 3s(^S) to the ground state (note that 
the ground state is split into three j sub-levels and appropri- 
ate weights have to be used; also the radiative rates of these 
j sub-state transitions differ) . The weights of the upper and 
lower states are (J'3s(f5s),2p(3p) ; fssC^S) is the normalised ef- 
fective recombination fraction into the 3s(^S) state and C 
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is the de-excitation rate due to thermal electrons which is 
much smaller than A/t. The 3s(^S) to 2p(^P) transition is 
spin-forbidden but dipole-allowed and therefore coUisional 
rates are w eaker than radiative r ates by several orders in 
magnitude (|van Regemorte il 11962') for typical SN nebular 
densities. For noC > 10^ cm~^ this assumption breaks down 
since C and r will increase with density and clumping factor. 

Equation[TT]can be used to calculate the Sobolev optical 
depth of the 7774 A line 

-^7774iS'3p(=P)^7774n3g(5s)[C] 
T-3.(55) ^ , n3,(5s) > n3p(5p) 



(14) 



with the transition wavelength A7774 = 7774 A, t the epoch 
in seconds, ^47774 the radiative rates from the 3p(^P) to 
the 3s(^S) state (note that the 3p(^P) state is split into 
three j sub-levels and appropriate weights have to be used) 
and ff3p(5p)^3s(5s) the weights of the upper and lower states. 
Again we compare this estimate with numerical results from 
our nebular code in Appendix A. 

At this point it is important review the concept of 
clumping in more detail. It is assumed that the ejecta are 
distributed into small "blobs" with size much smaller than 
the typical scale of the SN, covering the SN volume ho- 
mog eneously (e.g. iLi fc McCra^ll992l ). Therefore, although 
the optical depth is increased locally in the emitting region, 
the global optical depth of the SN remains constant. This 
means that, when scattering radiation from remote regions 
of the SN, the opacity is not influenced directly by clump- 
ing. On the other hand, when scattering radiation in the 
emission region, the optical depth increases proportionally 
to the clumping factor. 

In addition, there is an increase of the ratio ^'^''^^ 

' "Tot 

(see Equation llip . which results from the increased self- 
absorption of ground state transition radiation. Therefore, 
even when scattering radiation from remote regions, clump- 
ing increases the optical depth and so it does influence the 
7774 A resonance scattering of background radiation. 

The 3p(^P) state is separated from the 3s(^S) state by 
~ 1.6 eV, therefore thermal electron excitation is possible, 
but the population of the 3s(^S) state by thermal electron 
excitation is too low. Therefore the state must be populated 
by recombining electrons. 

The thermal electron excitation coefficient from 3s(^S) 
to 3p(^P) is given by 



1^ o c in-6 "eC ^7774 / -E7774 s 

C3p(5p) ~ 8.6 X 10 exp — 

rV2g3^(5gj kT (15) 



where all constants have the same meaning as in Equation 
[141 £7774 is the energy correspo nding to 7774 A and 07 774 
= 25.1/40.2 (at 5000/10000 K: iBhatia fc Kastned ll995'l is 
the efi'ective collision strength of the O I 7774 A transition. 
One can compare the population of the 3p(^P) state due 
to thermal electron coUisional excitation and recombination 



respectively by calculating their ratio 

C3p{5p)»^3s(5S) 



3p{5P) - 



yo[I]'^o[Il,f3p(Sp) 

(l + 7^)Co7 — ^ — ^("Tot - no[i])no[i]C 



(16) 



7; 



3p(5p) 



The ratio of the O I 7774 A luminosity induced by thermal 
electron excitation compared to the total deposited luminos- 
ity is then given by 

^3p(Sp) n -U-P^2Co^7774/3s£S)20 Upy^ 3 2 

(1 + K) 1 (1 - )'^o[iiC 



Wo[i]A 



nTot 



~ 10"^'f[days](l - ^^)r^^mC^ T = 5000 K 
riTot 

(17) 

This estimate becomes invalid for no[i]C > 10^ cm~'^ since 
Equation[TT]had been used for its derivation. Equation [T7] is 
compared to results from our nebular code in Appendix A. 
Whether coUisional excitation of the O I 7774 A line is 
negligible depends mainly on the temperature, which de- 
creases with time, on the density of neutral oxygen, which 
decreases with t~^, and on the clumping factor. Therefore, 
especially for early epochs (high temperature, high density) 
and for large clumping factors ((" ^ 1) the O I 7774 A 
line may be excited by thermal excitation of recombin- 
ing electrons. Around 150 days the core (oxygen) density 
of CC-SNe is typically of order 10^~* cm~^, which gives 

~ (lO"'^ to 10"^)C^ assuming a temperature of 5000 



K. 

From the above considerations it becomes clear that 
clumping influences the recombination line strength oifly 
weakly, but increases the effect of resonance scattering and 
thermal excitation of the O I 7774 A transition strongly. For 
both processes the 3s(^S) population is provided by recom- 
bination and not by thermal electron excitation from the 
ground level. 



2.4 Emission versus absorption line shapes 

Since the [O l] AA 6300,6363 doublet has often been used 
to probe the ejecta velocity and the geometry of CC-SNe 
cores (see Section [TJ it is interesting to study the profile of 
the [O l] 7774 A fine. If this line were caused by recombina- 
tion or thermal excitation alone, one would expect that the 
shapes of the 6300, 6363 and 7774 A lines should be at least 
approximately similar. 

However, as we have shown above, O I 7774 A might 
result from line scattering even at 200 days. At later times 
the line usually becomes very weak and it is difficult to use 
it to probe the ejecta geometry. The profile of a scattering 
dominated line depends on the fiux distribution, so one can- 
not expect that the shape of the 6300, 6363 A lines and the 
7774 A line will be even approximately similar. 

To demonstrate this, we use a multi-dimensional Mon- 
tecarlo code. A photon background is generated and 
resonance-scattered in an (asymmetric) oxygen distribution 
(e.g. taken from a nebular model; the nebular code in its 
current version cannot simulate scattering processes). We 
calculate the [O l] AA 6300,6363 doublet for an asymmetric 
oxygen model and use the 7774 A line opacity calculated in 
the nebular code to compute the absorption of a broad-band 
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Figure 4. Left side: flux of a synthetic [O l] AA 6300,6363 doublet 
for an asymmetric model (double-peaked shape). Right side: the 
7774 A line opacity calculated from this asymmetric model was 
used to calculate the flux of the absorption line at 7774 A (dashed 
vertical line). The flux is normalised to an arbitrary constant. 
The formation of the absorption line strongly depends on the 
assumed background (indicated by the dashed horizontal line; it 
not only depends on the absolute background flux but also on the 
background rest-frame wavelength which determines the emission 
site) and therefore the absorption line shown would look different 
for a different background radiation fleld. The asymmetry has 
some influence on the absorption line, but the clear double-peaked 
shape of the [O l] AA 6300,6363 doublet is not seen. 



background (3000 — 10000 A) emitted in the inner region 
of the model SN. The formation of the absorption line is 
very sensitive to this background and therefore our results 
(see Figure[4]) illustrate just one of many possibilities. As ex- 
pected, the absorption line has quite a different shape than 
the emission line. 



2.5 lonisation Rates for the nebular code 

To improve the accuracy of the recombination line calcu- 
lation in our nebular code we updated the treatment of 
ionisation for all the ions treated. This is necessary, since 
the formation of recombination lines is much more sensitive 
to the ionisation rate than the formation of other emission 
lines, and higher accuracy is needed. In the older version of 
the nebular code the ionisation rates were calculated using 
a simple analytical form. 

To obtain more accurate ioni sation rates we calculate 
the average 'work per ion ' W (see lAxelrodlll98dl , which is 
the amount of energy lost by a non-thermal electron (with 
energy 3> 1 keV) to cause one ionisation. The ionisation rate 
is then computed dividing the total deposited luminosity by 
this energy. 

The value of W can be calculated averaging the 'work 
per ion' W* of all individual ions over the total ionisation 
cross sections. W* can be obtained for any ion integrating 
the ratio of the energy-dependent ionisation cross section (Ji 



and the energy loss functions Li (|Axelrodlll98Cl l 

E 

max 



dE' 



(18) 



Highly energetic electrons lose most of their energy to 
ionisation and to the secondary electrons produced in the 
ionisation, but also by non-thermal excitation. At low ener- 
gies, losses to the electron plasma become important, how- 
ever most of the primary energy is lost in atomic collision 
processes. 

lonisation cross sections and loss functions can be ob- 
tain ed for any ion from the literature at low electron energies 
(e.g. lLotj[l967l . ll970l : lLennon et al.lll988l '). At high electron 
energies they can be compute d by appropri ate extrapola- 
tions of the low energy values (|Axelrodlll980l ). 

The ionisation rate of any ion is then given by 

LocpO'i _ i/Dcp 



Yi 



(19) 



NTotW <a> iVTotWi 

The ionisation rates mainly affect the recombination 
lines but also have some influence on all other emission lines. 
This will be discussed in Section [S] 



3 A SHELL MODEL OF SN 2002AP 

SN 2002ap is classified as a broad-lined SN of Type Ic. It 
had an ejecte d mass of ~ 2.5 a nd a kinetic energy of ~ 

4 • lO'^^ ergs (|MazzaH et all 12007^ 1. The distance and red- 
dening to SN 2002ap are only known approximately. To be 
consistent with previous work we us e ^ = 29.50 mag an d 
E(B - V) = 0.09 mag as done by iMazzali et al.1 l|2002h : 
lYoshii et al] (|2003l '): iMazzali et"al] (|2007al 'l. It is not clear 
whether SN 2002ap was a sph erical symmetric event or not. 
Recently. iMaurer et al.l l|2009l ) have shown that an asymme- 
try might be observable in all broad-lined CC-SNe, however 
this point is not clear yet. Th e spectra of SN 2002a p used 
in t his work are those use d bv IMazzali et al] l|2007al ). 

IMazzali et al.l (|2007al ) found appropriate one-zone and 
shell models of the SN 2002ap nebular ejecta. We quickly 
summarise the main findings for the shell models here. Using 
a filling factor of 0.1 (C = 10) in the ^^Ni rich regions, a total 
ejecta mass of ~ 2.5 M© was found, containing roughly 0.11 
M0 ^"^Ni and 1.3 M© of oxygen. 

In contrast to this previous work, which aimed at mod- 
elling each spectrum individually, here we try to find one 
smg le ^^m/O model which can produce the time evolution 
of the different spectra at all observed epochs consistently. 
Special attention is paid to reproducing the exact shape of 
the [O l] AA 6300,6363 line profile, which is a tracer of the 
distribution of oxygen, the most abundant element of the 
nebula. 

Using th e same distance modulus and reddening as 
IMazzali et al.l (|2007al ) and clumping factors of 5 and 25 [we 
had to use a clumping factor of 1 for the innermost shell 
of the = 25 model to avoid the fo rmation of sharp high 
density lines which are not observed; IMazzali et al.l (|2007al ) 
had used a clumping factor of 10; this will be discussed in 
Section [5] we obtain models similar to the ones given in 
IMazzali et al.i ((2007al ) giving reasonable agreement with the 
observations at all epochs (see Figures [5] and [7|. The forbid- 
den oxygen lines are reproduced at all epochs by one and 
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the same oxygen distribution, which provides evidence for 
the reliability of this model. 

Our '^''Ni zone with a total mass of 1.2 M0 extends 
out to 800 km/s, containinR 'y 0.07 M© of ^"^Ni (in agree- 
ment with lMazzaU et al.|[20o3 l and ~ 0.7 M© of O. More 
mass is located at higher velocities, but the nebular mod- 
elling becomes inaccurate for the outer regions. The exact 
values depend on the clumping factor as well as on the ^''Ni 
distribution, which is not known. 

The O I 7774 A line is excited by recombination and 
thermal electron excitation in our nebular modelling (no line 
scattering) and is too weak to explain the observations at 129 
and 163 days using low clumping factors ((" ~ 5). The ^ = 25 
model can reproduce the observations of the O I 7774 A line 
at 129 and 163 days better than the = 5 model. This 
results from thermal electron scattering of 3s(''S) electrons, 
which are provided by recombination. In general the C = 5 
model reproduces the formation of the forbidden lines better 
than the = 25 model at early epochs, since high density 
lines, which are not observed, form in the = 25 model. 

At days 192 and 229 both models produce too little flux 
at 7774 A, which is probably influenced by line scattering 
at these epochs (the optical depth of the inner shells of the 
O I 7774 A line is high for both models at these epochs). 
The scattering process is not simulated and therefore the 
synthetic flux is too low. 

At later epochs the O l 7774 A line is consistently repro- 
duced by recombination alone in both models when taking 
into account the uncertainty due to the background around 
7774 A, which is not reproduced by the nebular code. 

The optical depth of the 7774 A line is shown for all 
inner shells (five shells between and 5000 km/s) of our 
(" = 5 model in Figure [6] It does not fall below one before 
210 days, indicating that there can be strong line scattering 
up to ~ 200 days depending on clumping (for the C, = 2^ 
model the optical depth is higher and the O I 7774 A line 
becomes optically thin later). 

For clumping factors below 5 neither line scattering nor 
thermal excitation was strong enough to allow sufficient for- 
mation of the O I 7774 A line before ~ 200 days. 



4 A 2D MODEL OF SN 1998BW 



SN 1998bw is classified as a broad-lined SN of type Ic 
(ejected mass ~ 14 and kinetic energy ~ 6 ■ 10^^ ergs; 
iNakamura et all I2OO0I '). We use a distance modulus of /i 
= 32.89 and an e xtinction of Av = 0.2 mag as given by 
iPatat et al.l l|200ll ). SN 1998bw was accompanied by the low 
energy, long-duration GRB 980425. Due to this GRB-SN 
connection there are speculations that SN 1998bw was a 
highly aspherical event. The se speculations are supported by 
polar i sation measureme nts (|Kav et al.ll 19981 : llwamoto et al.l 
Il998l : iPatat et aklbOOlll. which where interpreted as a SN 
axis ratio of 2:1 jHoflich et all 1 19991 ). Some indication for 
core ejecta asymmetry ha d also been found i n the nebu- 
lar spectra of SN 1998bw (|Mazzah et al.ll200ll ). The spec- 
tra of_SN_J^98bw used in this work are the s ame as used 
by Mazzali et al.] (|200ll ) originally presented bv I Pat at et al.l 
l|200ll ). 

The nebular phase of SN 1998bw has been well stud- 
ied by means of nebular modelling l|Mazzali et al.l I2OOII : 




200 250 500 
Epoch [d] 



Figure 5. SN 2002ap, models for a clumping factor = 5. Left 
side: normalised nebular spectra (from top to bottom: 129, 163, 
192, 229, 245, 253, 281, 343 and 394 days after explosion) are 
shown in black with synthetic spectra on top in red. Right side: 
the region between 7600 A and 8200 A is enlarged to enable a com- 
parison of observed (black) and synthetic (red) 7774 A flux. At 
epochs of 129 to 229 days, the simulated flux is too low to match 
the observations. Beginning at 245 days, considering a constant 
off set caused by background-flux which is not reproduced by the 
synthetic spectrum, the synthetic 7774 A line starts to become 
consistent with the observations. 




150 



200 



250 300 
Epoch [d] 



350 



Figure 6. SN 2002ap, models for a clumping factor f = 5. The 
logarithm of the 7774 A line optical depth in our model for the in- 
ner shells (from top to bottom: - 1000 km/s, 1000 - 2000 km/s, 
3000 - 4000 km/s, 4000 - 5000 km/s and 2000 - 3000 km/s). 
The 2000 — 3000 km/s shells contains very little oxygen in order 
to reproduce the narrow peak on top of the [O l] AA 6300,6363 
doublet lines and therefore it has low optical depth. Apart from 
the innermost shell, the optical depth drops below one around 
220 days. This optical depth is calculated for scattering of re- 
mote emission radiation, which has a weaker dependence on the 
clumping factor than for local scattering. 
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Figure 7. SN 2002ap, models for a clumping factor = 25. Left 
side: normalised nebular spectra (from top to bottom: 129, 163, 
192, 229, 245, 253, 281, 343 and 394 days after explosion) are 
shown in black with synthetic spectra in red. Right side: the region 
between 7600 A and 8200 A is enlarged to enable a comparison 
of observed (black) and synthetic (red) 7774 A flux. At epochs of 
129 and 163 days, the simulated flux matches the observations at 
7774 A much better than for low clumping factors but in general 
the fit to the forbidden line emission is worse (e.g. the [O i] 5577 A 
line and the Ca II IR-triplet line at ~ 8500 A become too strong 
at early epochs). Around 200 days line scattering seems necessary 
to explain the observed 7774 A flux. 

iMaeda et al] 120061 ). While iMazzali et al.l |200j) attempted 
to model individual nebular spectra by one-zone models 
(concluding that this is insufficient), iMaeda et al.l (|2006l ) 
computed two-dimensional synthetic nebular spectra for 
models obtained from hydrodynamical simulations, finding 
good agreement with spectra of individual epochs. However 
no model has so far described the time sequence of nebular 
spectra of SN 1998bw consistently. 

In this section we try to find a single two-dimensional 
model which is consistent with the full time evolution from 
108 to 388 days after explosion. While successful models 
for individual epochs can be found even in one dimensional 
modelling, it seems impossible to model the time sequence of 
all spectra with a single one-dimensional model (in contrast 
to SN 2002ap, where this approach works well). 

Since SN 1998bw is quite massive compared to other 
CC-SNe (like SN 2002ap for example), the transition to 
the nebular phase occurs at rather late times, although the 
ejecta expand at high velocities. The spectra at 108 and 139 
days after explosion are not strictly nebular and therefore 
large deviations between the model and the observations are 
to be expected when trying to reproduce these spectra (es- 
pecially at 108 days) with our nebular code. 

When considering two dimensions one is immediately 
confronted with the problem that the parameter space is 
much larger [angular distribution of elements, observer in- 
clination]. In addition, the computation time to obtain one 
synthetic spectrum increases dramatically. Since hundreds of 
spectra have to be computed to obtain a good model, this 
technical problem can make modelling unfeasible. Therefore 
some assumptions about the SN geometry have to be intro- 
duced. 



There are several good arguments to assume that SN 
1998bw might consist of some kind of two-dime nsional 'jet -|- 
disc' s tructure, as described by for example bv lMaeda et al.l 
(|2006l) (see Section [1]). Furthermore, since the SN was ac- 
companied by a CRB, it seems likely that the observer incli- 
nation is not too far from polar, although this is uncertain 
since the opening angle of GRB 980425 is not known. 

Therefore we work with a parametrised two-dimensional 
model which consists of a polar zone (0° to 45°) and an equa- 
torial zone (45° to 90°). This simplification will introduce 
some error, since most likely it does not represent physi- 
cal reality. However, it seems sufficiently exact to obtain an 
acceptable fit at all epochs with a sing le 5^Ni/0 model. A 
model with more degrees of freedom is highly degenerate 
anyway, since the information that can be extracted from a 
series of spectra F(i',t) is limited. 

Because of this simplification the reliability of our model 
is in question, but this kind of uncertainty is inherent to all 
multi-dimensional modelling. We performed tests for view- 
ing angles of 0°, 15° and 30°. Using the smaller viewing 
angles the iron lines which are mainly produced in the jet- 
like structure become very narrow, which is not observed. 
Therefore we decided to use a viewing angle of 30°. How- 
ever this might be a consequence of our simplified geometry 
and therefore it is not possible to obtain valuable informa- 
tion about the observers inclination from our modelling ap- 
proach. 

Our model is r oughly consistent with previous findings 
l|Maeda et al.ll2006l ). The polar zone contains ~ 0.24 M© of 
®^Ni at velocities below 12000 km/s. The equatorial zone 
also contains ~ 0.24 Mq of ®^Ni (note that the equatorial 
zone has more than twice the surface area of the polar zone) , 
but located at velocities below ~ 8000 km/s. The total mass 
below 12000 km/s is estimated to be ~ 2.7 Mq containing a 
total oxygen mass (below 12000 km/s) of ~ 2 Mq. There is a 
lot of material at higher velocities, but the nebular modelling 
becomes inaccurate for the outer regions because the density 
is too low. 

Again we find that the O I 7774 A line is not repro- 
duced as due to rec ombination by our m odel with clumping 
C = 1 and C = 5. iMazzah et~al] l|200ll ) and iMaeda et all 
( 2006. ) had used a clumping factor of 10. This will be dis- 
cussed in Section[S] for the spectra at 108, 139 and 214 days 
(see Figures |S] and llOjl . At 349 and 388 days we can repro- 
duce the O I 7774 A line consistently with the observations. 
There is some background around 7774 A which is not re- 
produced by the nebular code sufficiently exactly. Shifting 
the oxygen recombination line upwards by the possible back- 
ground, the agreement seems reasonable for both clumping 
factors. However the large background-line ratio makes an 
exact comparison impossible. 

We tried to explain the formation of O I 7774 A at the 
earlier epochs (108, 139 days) by increasing the clumping 
factor, which in turn increases the thermal electron excita- 
tion from the 3s(^S) to the 3p(^P) level. However, this does 
not seem to work. Both [O i] 5577 A and the Ca ii IR-triplet 
become very strong at clumping factors > 10, while the syn- 
thetic flux around 7774 A increases but seems still too weak 
to explain the observations. Owing to the high density of 
SN 1998bw the nebular approach might be not suitable for 
such early epochs (the strong continuum which is observed 
is not reproduced). In addition there are uncertainties on 



© 0000 RAS, MNRAS 000, 000-000 



10 Maurer et al. 




4000 5000 6000 7000 8000 9000 7600 8200 
Wavelength [A] 




200 250 500 
Epoch [d] 



Figure 8. SN 1998bw, case of no clumping {C, = 1). Left side: nor- 
malised nebular spectra of SN 1998bw (from top to bottom: 108, 
139, 214, 349 and 388 days after explosion) are shown in black 
with the synthetic spectra in red. On the right side the region 
between 7600 A and 8200 A is enlarged to enable a comparison of 
observed (black) and synthetic (red) 7774 A flux. The two earliest 
spectra (108, 139 days) of SN 1999bw are clearly not nebular and 
the agreement between observations and simulations is expected 
to be poor. However, the evolution of total and oxygen 6300 A 
flux as well as the shape of the [O l] AA 6300,6363 doublet line are 
reproduced at all epochs much better than by a one-dimensional 
model. At epochs of 108 to 349 days, the simulated flux is too 
low to match the observations of the 7774 A line. At 388 days, 
considering the constant offset caused by some background flux 
which is not reproduced by the synthetic spectrum, the synthetic 
7774 A line seems to be consistent with the observations. Unfor- 
tunately the noise level is high and a detailed comparison is not 
possible. 



the atomic data (especially the collision strengths), which 
may influence our conclusion. 

In Figures [91 and [TT] we show the opacity of the 7774 A 
line for the flve innermost shells for clumping factors of ^ = 
1 and ^ = 5. There will be 7774 A line absorption up to 
~ 200 — 300 days depending on clumping. 

Again, at least moderate clumping ((" ~ 5) seems nec- 
essary in order to provide enough O I 7774 A line opacity to 
reproduce the observations at 214 days. 



5 DISCUSSION 



In Section [2. II we have obtained normalised effective recom- 
bination fractions using available atomic data and comput- 
ing missing atomic data. The agreement with atomic data 
available in the literature is go od. Effective re c ombin ation 
rates had been given before by ljulienne et al.l (|l974h at a 
temperature of 1160 K. These results are reproduced well. 
We extend the temperature range up to 10000 K, which 
should be an upper limit for temperatures of SN nebular 
ejecta between 100 and 600 days. Typical uncertainties in 
the atomic data seem to have small influence. 

In Sections 12.21 and 12.31 we have obtained estimates for 
the formation of oxygen recombination lines and line absorp- 
tion by the 3s(^S) state of neutral oxygen. We have shown 
that the luminosity of recombination lines is weak at any 



Figure 9. SN 1998bw, case of no clumping (f = 1). The logarithm 
of the 7774 A line optical depth in our model for the inner shells 
(from top to bottom: - 1000 km/s, 1000 - 2000 km/s, 2000 
- 3000 km/s, 3000 - 4000 km/s and 4000 - 5000 km/s). Solid 
lines show the equatorial region, dotted lines the polar region. 
Apart from the innermost shell, the optical depth drops below 
one at ~ 210 days. This optical depth is calculated for scattering 
remote emission radiation, which has a weaker dependence on the 
clumping factor than for local scattering. 
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Figure 10. SN 1998bw, case of clumping {C, = 5). Left side: 
normalised nebular spectra (from top to bottom: 108, 139, 214, 
349 and 388 days after explosion) are shown in black with the 
synthetic spectra in red. Right side: the region between 7600 A 
and 8200 A is enlarged to enable a comparison of observed (black) 
and synthetic (red) 7774 A flux. The two earliest spectra (108, 139 
days) of SN 1999bw are clearly not nebular and the agreement be- 
tween observations and simulations is expected to be poor. How- 
ever, the evolution of total and oxygen 6300 A flux as well as 
the shape of the [O i] AA 6300,6363 doublet line are reproduced 
at all epochs much better than by a one-dimensional model. At 
epochs of 108 to 349 days, the simulated flux is too low to match 
the observations of the 7774 A line. At 388 days, considering the 
constant offset caused by some background flux which is not re- 
produced by the synthetic spectrum, the synthetic 7774 A line 
seems to be consistent with the observations. Unfortunately the 
noise level is rather high and a detailed comparison is not possible. 
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Figure 11. SN 1998bw, case of clumping (C = 5)- The logaritlim 
of tlic 7774 A line optical depth of our ejecta model for the inner 
shells (from top to bottom: - 1000 km/s, 1000 - 2000 km/s, 
2000 - 3000 km/s, 3000 - 4000 km/s and 4000 - 5000 km/s). Solid 
lines show the equatorial region, dotted lines the polar region. 
Apart from the innermost shell, the optical depth drops below 
one at ~ 310 days. This optical depth is calculated for scattering 
remote emission radiation, which has a weaker dependence on the 
clumping factor than for local scattering. 



epoch. The influence of clumping on the recombination line 
is weak as well. Clumping will increase the recombination 
rate, but will decrease the number of O ll ions at the same 
time. The total number of recombinations, which is equiv- 
alent to the total number of ionisations, will increase only 
slightly. Therefore clumping offers no direct way for increas- 
ing O I 7774 A emission considerably. 

We have shown that there is line scattering of 7774 A 
photons even at very late epochs. This is possible because 
of the large number of ground state O I atoms, which res- 
onantly scatter the UV radiation emerging from 3s to 2p 
transitions and slow down the de-population of the 3s state. 
Combined with recombination this leads to a high popula- 
tion of the 3s level, which in turn can resonantly scatter 
7774 A radiation. We have shown that, in contrast to the 
direct recombination line, the scattering line is very sensitive 
to clumping of the ejecta. 

In addition, for high clumping factors, there may be 
thermal electron coUisional excitation of the O I 7774 A line 
at early epochs. However, it is in question whether these high 
clumping factors are realistic since high density lines, which 
are not observed, may form. Uncertainties in the atomic 
data, especially the collision strengths, may influence our 
results. 

The considerations regarding the O I 7774 A line are 
also in principle valid for the [O l] 8446 A line [3s(^S) to 
2p(^P)] with the difference that the effective recombination 
rate into the 3s triplet state is smaller than into the quin- 
tet state. The optical depth of the 8446 A line is therefore 
expected to be lower, which is still sufficient to cause line 
scattering at very early epochs. However, the 8446 A line is 
blended with the blue wing of the [Ca ll] IR-triplet, which 
can explain why no [O l] 8446 A line is identified in SNe 
1998bw and 2002ap. 

We have shown that there is likely to be no other effec- 



tive excitation mechanism than 'recombination -I- thermal 
excitation' or 'recombination -f absorption' for O I n 3 
levels. Thermal electron excitation of n ^ 3 levels from the 
ground level is ineffective at typical nebular temperatures. 
Non-thermal electron excitation can also be ruled out. First, 
it is too weak and second, it would produce a temporal be- 
haviour of the O I 7774 A line which is not observed. As de- 
tailed modelling showed (Sections[3]and[4l) absorption might 
be too weak to explain the observations as long as there is no 
clumping of the ejecta ~ 5). To enable thermal electron 
excitation even larger clumping factors seem necessary. 

The concept of clumping has ofte n been used in stel- 
lar wind and SN rad iation physic s (e.g. lLi fc McCrav"l99' 
1993; Maz zali et alJ 2001; Macda etal.l 12006 : Mazzafi ct al. 
[2007a) without quantitative understanding of the physics 
behind eject a clumping, therefor e playing the role of a fit- 
ting factor. iLi fc McCravl (|l992l ') have modelled the ratio 
of the [O i] AA 6300, 6363 doublet lines and iLi fc McCravl 
(1993) the Ca ll emission of SN 19 8 7A using a clumpin g 



factor of ^ 1 . iMazzali et all (|200ll '): iMaeda et al] l|2006f) : 



iMazzali et al.l (12007^ ) have modelled the nebular spectra of 
SNe 1998bw and 2002ap finding that better fits to obser- 
vations are obtained using clumping factors of ~ 10 than 
when using no clumping. Our results are similar, but show 
some difference. We obtain good fits to all forbidden lines 
using clumping factors of 1 and 5 as well. The abundance 
of elements will change slightly depending on the value cho- 
sen. We even found a slight decrease in the quality of the fit 
when going from = 5 to ^ = 10 and therefore we used the 
lower value. This difference probably results mainly from 
our more accurate treatment of ionisation which increases 
the ionisation rate of carbon and oxygen considerably. This 
means that the electron density is increased, which can be 
mimiced by increasing the clumping factor and using lower 
ionisation rates. It might also be partially attributed to the 
fact that there is no strict criterion for the quality of a fit 
for nebular spectra. 

Still, the constraints we derived for the clumping factor 
of SNe 1998bw and 2002ap from the O I 7774 A scatter- 
ing line compare well to previous findings. We found that 
a clumping factor of ~ 5 seems necessary to provide suffi- 
cient opacity at 7774 A. Clumping factors of > 10 make it 
possible to model the O I 7774 A emission at early epochs 
without line scattering, but may cause the formation of high 
density emission lines at early epochs (100 — 150 days). One 
example is the [O l] 5577 A singlet line, which is observed 
to be weak. Since there are uncertainties on the atomic data 
it is not clear whether thermal electron excitation is impor- 
tant and whether the O I 7774 A line at early epochs is 
dominated by line scattering or thermal electron excitation 
of recombining electrons. 

At later epochs, between 150 and 250 days, the line is 
most likely a combination of line scattering and recombi- 
nation radiation, since thermal electron excitation becomes 
too weak for any reasonable value of clumping. Around ~ 
250 days the O I 7774 A line seems to become a true re- 
combination line where the fiux is mainly provided by the 
electrons cascading into the 3p(^P) level. 

In Section 12.41 we have studied the expected profile of 
the [O l] 7774 A line in light of our previous findings. As ex- 
pected, the absorption and emission line shapes caused by a 
certain ejecta distribution can show important differences. 
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Further, the profile of the absorption line depends on the 
background radiation field and one should be careful when 
using the shape of a scattering line for any kind of argumen- 
tation about ejecta geometry as long as it not clear how this 
line is formed in detail. 

In Sections [3] and |4] we obtained core ejecta models for 
SNe 1998bw and 2002ap. The primary goal of this modelling 
was not to re-derive ejecta properties, but to obtain reliable 
models for studying oxygen recombination consistently with 
other line formation. 

As a by-product we obtained estimates for the ^^Ni and 
oxygen masses of the cores of these SNe. We find that SN 
2002ap can be described very well by a one-dimensional shell 
model (as already found before), which does not necessarily 
mean that there is no asymmetry. Our ^^Ni (~ 0.07 M©) 
and oxygen mass estimates seem consistent with previous 
work. For SN 1998bw one-dimensional modelling seems to 
be insufficient to obtain an acceptable fit at all epochs using 
a single model. An acceptable fit is obtained with a two- 
dimensional 'jet -|- disc' model (this had been found before) . 
The fit could certainly be improved by increasing the degree 
of freedom of the model, but this would become extremely 
time-consuming. The model found in this work reproduces 
the observations at all epochs with increasing accuracy at 
later times, which is expected. The estimate of the ^^Ni (~ 
0.48 M0) and oxygen mass are consistent with previous re- 
sults. 

Using these ejecta distributions we were able to repro- 
duce the formation of O I 7774 A at late epochs. There is 
enough line opacity at 7774 A at early and intermediate 
nebular epochs to produce a strong scattering line. Since 
the background radiation field is important for the calcula- 
tion of this absorption line, we were not able to show that 
our ejecta distribution reproduces the observations at early 
epochs exactly, however this seems likely. If there is strong 
clumping the O I 7774 A line may be additionally excited 
by thermal electrons. The oxygen distribution inferred from 
forbidden line observations allows the reproduction of the 
allowed oxygen lines at late epochs. Since different physics 
are involved, this provides a test for the consistency of the 
nebular modelling approach, which is passed for SNe 1998bw 
and 2002ap. 

Several other ions may produce absorption lines at 
wavelengths between 4000 and 10000 A (the part of the spec- 
trum that is usually observed). Since the abundance of these 
elements is much lower than that of oxygen, there will be no 
excited states with sufficient population to cause significant 
line scattering (an exception might be the Ca 11 IR-triplet). 
Allowed ground state transitions with energies between ~ 1 
and 3 eV are interesting since their optical depth can be high 
owing to the high radiative rates of allowed lines and to the 
possibly sufficiently high ground state populations of these 
low abundance elements. The elements of interest are Na I, 
Mg I and Ca 11. These three elements have low quantum-level 
transitions with wavelengths of ^ 5890 A (Na l) , 4570 A 
(Mg l) and 3950 A , 7300 A and 8500 A (Ca 11), which can 
be present in sufficient amounts in SN nebulae. Iron group 
element lines, especially Fe 11, may be optically thick as well. 
However, since these ions are complex we do not treat them 
explicitly here. 

If 0.01 M0 of Na I were to be distributed homogeneously 
within 5000 km/s, the optical depth of the Na I 5890 A 



line would be ~ 6 x 10^ at 200 days (this result is obtained 
by balancing thermal-electron excitation with radiative and 
coUisional de-excitation analogously to Section [2. 2|l . scaling 
linearly with Na I mass and with epoch . 

If 0.01 Mq of Mg I were to be distributed homoge- 
neously within 5000 km/s, the optical depth of the Mg l] 
4570 A line would be ~ 1 at 200 days, scaling linearly with 
Mg I mass and with epoch t~^. 

Three Ca 11 lines might be seen in absorption even at 
late epochs (~ 3950 A, 7300 A and the IR-triplet 8500 A). 
To estimate their strength, it is important to know the rel- 
ative population of the Ca 11 3d to 4s states, which can be 
approximated by 



'^4s ^3d_»4s + C'3d->4s 

where C4s->3d and C3d->4s are the coUisional rates for the 4s 
to 3d transition and A' is the radiative rate from 3d to 4s, 
reduced due to resonance scattering. 

If 0.01 Mq of Ca II were to be distributed homoge- 
neously within 5000 km/s, the optical depth of the Ca 11 
3950 A and [Ca 11] 7300 A lines would be ~ 2 x 10^ and 
~ 0.02 respectively, at 200 days, scaling linearly with the 
Ca II mass and with epoch t~^. The optical depth of the 
Ca II IR-triplet depends on the population of the 3d state. 
Using the 7300 A optical depth estimated above, the optical 
depth of the Ca 11 8500 A transition is ~ 430 (rie = lO'' and 
T = 5000 K) at 200 days, with a stronger dependence on 
Ca II mass and epoch than the opti cal depths at 3950 A and 
8500 A (see also lLi fc McCravlll993l . for a detailed discussion 
of Ca II line formation) . 

The exact mass of Na I, Mg I and Ca 11 is not well 
known, since estimates of their masses are not very accurate. 
However, typical masses in CC-SNe cores for Na, Mg and Ca 
might be of the order of 0.01 M©. 

Therefore, even at late times (> 200 days) there is 
strong line scattering at wavelengths between ^ 4000 and 
10000 A. The Ca 11 3950 A, Na I 5890 A and O l 7774 A lines 
cause strong scattering in any reasonable scenario. The [Mg 
l] 4570 A, O I 8446 A and Ca 11 8500 A lines may have suffi- 
cient optical depth for line scattering, depending on epoch, 
degree of ionisation, total mass and clumping. 



6 SUMMARY AND CONCLUSION 

We have computed temperature-dependent effective recom- 
bination rates for neutral oxygen in a temperature range 
suitable for all types of SN nebulae (at epochs between 100 
and 600 days). Since oxygen is the most abundant element 
in stripped-envelope CC-SNe, oxygen lines are of special in- 
terest among other recombination lines. 

We obtained core ejecta models for CC-SNe 1998bw and 
2002ap. Similar models had been derived previously. Using 
these oxygen profiles, the O I 7774 A recombination line, 
which is the strongest observed recombination line, is calcu- 
lated and compared to observations. We show that up to late 
epochs pure recombination is too weak to power O I 7774 A. 
At earlier epochs the line is powered by scattering and possi- 
bly by thermal electron excitation of recombining electrons. 
In both scenarios the population of the 3s(^S) state by re- 
combining electrons is the key to the O I 7774 A emission. 
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We derived estimates for the strength of the O I 7774 A 
hric resulting from oxygon recombination, for the time- 
dependent optical depth of this line and for excitation of re- 
combining electrons by thermal collisions. These estimates 
give insight into the formation of the O I 7774 A line in 
CC-SNe. We have shown that while the recombination line 
strength depends on clumping very weakly, clumping does 
have a strong influence on the absorption and thermal elec- 
tron excitation of this line. 

Constraints on the clumping factor have been rare so 
far. A common choice is a value of ^ ~ 10, consistent with 
the findings of this paper. At the current level of accuracy it 
is only possible to set upper and lower limits to the clumping 
factor (100 > C > !)• 

Our results imply that the [O l] 7774 A line should not 
be used as a tracer of the core ejecta before 250 days, unless 
one explicitly models the background radiation field and the 
recombination and absorption processes involved. Since the 
clumping factor can not be obtained with high accuracy and 
the background around 7774 A seems to consist of a super- 
position of several weak lines this could become a difficult 
task. It is not clear which elements or ions produce the ob- 
served background flux in this region. Iron group elements 
seem to be promising candidates, but currently the accuracy 
of our calculations is not sufficient. 
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Figure 12. The ratio 4y^[i]i?o[i]»^TotC/(-Ro[l]"FeC + 5^0[l])^ for 
different clumping factors (C = 1, 10, 20, 30, 40, 50, 60, 70, 80, 
90 and 100 from top to bottom) at 1000 days calculated for the 
same model as used in Figures I3I13I and 1141 For this model the 
ratio is much larger than one for any reasonable clumping factor 
(1 - 100) at epochs earlier than 600 days. 
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APPENDIX A 

In Section 12.21 and 12.31 we have derived several estimates 
for quantities relevant for calculating O I 7774 A nebular 
emission. To derive these estimates, we used some approx- 
imations. These are justified here by comparison to results 
from the nebular code, which are not based on these approx- 
imations and will show deviations as soon as our approxi- 
mations fail. An important difference is introduced by the 
temperature. While the electron temperature is calculated 
considering hundreds of emission lines in the nebular code, 
this is not possible for our estimates and we use a tempera- 
ture of 4000 K for the calculations. 

To perform these tests we set up a test one-zone model 
with a total mass Miot ~ 0.5 Mq, an expansion velocity of 
5000 km/s, containing 80% oxygen, 10% ^^Ni, ~ 9% carbon, 
and other elements like calcium, magnesium and sodium. 

Usually one can assume 4Fo[j]7?o[i]'^TotC/(-Ro[i]'^FcC + 
^o[i])^ ^ 1, which depends on the iron to oxygen ratio, 
the deposited luminosity, the electron temperature and the 
clumping factor. In Figure[T2]we show this ratio for our test 
model using clumping factors of one to one hundred. At any 



Figure 13. The solid line is a computation of Equation |9] per- 
formed for the O I 7774 A line assuming a constant temperature 
of 4000 K at all epochs for the one-zone model described in the 
text. The luminosity of this line is compared to the total deposited 
luminosity. The results are given in %. The recycling fraction was 
set to the ratio of ground state to total recombination rate at 
4000 K. The crosses show data points obtained with the nebular 
code for the same model, where the temperature is computed at 
all epochs self-consistently and ionisation is treated in detail. Con- 
sidering the large change of the temperature between 100 and 600 
days (6800 K and 2800 K respectively) the agreement seems rea- 
sonable (note that the temperature calculated by the code drops 
below 4000 K around 450 days). The minimum around 200 days 
is present in both curves. 
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Figure 14. The solid line shows a computation of Equation 1141 
performed for the logarithmic line opacity of the O I 7774 A line 
assuming a constant temperature of 4000 K at all epochs for 
a one-zone model described in the text. The crosses show data 
points obtained with the nebular code for the same model, where 
the temperature is computed at all epochs self-consistently and 
ionisation is treated in detail. The agreement is good, especially 
around 450 days, when the simulated temperature drops below 
4000 K, which we assumed in our calculation. The line optical 
depth is greater than one until day ~ 200. 
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Figure 15. The O I 7774 A line luminosity from pure recombina- 
tion and thermal excitation of recombining electrons relative to 
the total deposited luminosity in % for the test models described 
in the text using clumping factors between 1 and 100. Solid lines 
show estimates obtained from Equat ions l9land ll7l for oxygen den- 
sities of lO'', lO'^'S and 10* cm from bottom to top. Results 
from our nebular code are shown by diamonds (10^ cm~^), tri- 
angles ilO"^-^ cm~^) and squares (10* cm~^). The agreement is 
good. Small deviations are predominately caused by our rough ap- 
proximation of the temperature. Filled squares mark rioC ~ 10' 
cm~'^ where our estimates are expected to become inaccurate. At 
noC > 10^'^ cm~^ our approximation breaks down completely, 
but this regime is probably not important for the nebular phase 
of stripped CC-SNe. The 7774 A luminosity increases with the 
increasing clumping until the density becomes high enough to 
depopulate the 3s(^S) state by thermal electron collisions effec- 
tively. 

epoch of interest (100 to 600 days) and for any clumping 
factor the assumption is valid. The assumption may fail for 
clumping factors ~ 1 at lower densities and for higher ^®Ni 
to oxygen ratios. For clumping factors C ^ 1 the assumption 
will hold for any plausible CC-SN scenario. 

In Figures [TH] and [T?] we compare Equations [5] and [TJ] 
(temperature set to T = 4000 K) with the results obtained 
from our nebular code using our test model. The tempera- 
ture is kept constant in the calculation, while it is calculated 
self-consistently in the nebular code, influencing the recom- 
bination fraction and recombination rate. Also the ionisation 
of ions other than O I and Fe I has been neglected in the cal- 
culation. Despite the approximations, the agreement seems 
reasonable. 

In Figure [15] we compare the O I 7774 A flux from pure 
recombination and electron scattering of recombining elec- 
trons (Equations [5] and I17[) to results obtained from our 
nebular code using variations of our test model. We set the 
epoch to 150 days and use oxygen densities of 10^, 10^'^ 
and 10* cm~'^, which should cover the range expected for 
stripped CC-SNe. The oxygen fraction is again set to 80%. 
For our estimates we use a temperature of 5000 K (the 
temperature calculated by the nebular code varies between 
4500 K and 6300 K for these models). The breakdown point 
of our estimate (noC > 10^ cm~^) is indicated by filled 
squares. Up to the breakdown point the agreement is good. 
For larger clumping factors the agreement is acceptable for 
noC < 10^'^ cm""'. For larger clumping factors our estimate 



over predicts the 7774 A luminosity. In this regime the popu- 
lation of the 3s(*S) level is no longer controlled by the effec- 
tive radiative rate but rather by collisional de-excitation by 
thermal electrons. It seems unlikely that such high density 
regimes can play a role in CC-SN nebulae. 



APPENDIX B 

Normalised effective recombination fractions for triplet and 
quinted states of O I n = 3. For the triplet states we show 
calculations for the cases of optical thin * and thick ground 
states. These numbers have to be multiplied by the total ra- 
diative recombination rate to obtain the effective recombina- 
tion rates. Since the total radiative recombination rate also 
contains a direct ground state component and since lower 
levels contain cascade contributions from the higher levels 
these fractions do not add up to one. 

At high densities, collisional transitions become impor- 
tant, which are not included in our cascade calculation. At 
densities We < 10^" cm~^ and temperatures of T ~ 5000K 
the error on the recombination rates owing to this effect 
should be less than 10%. 

At high temperatures di-electronic recombination be- 
comes important, while at low temperature and high densi- 
ties collisional recombination can become important. Under 
nebular conditions both effects should be of the order of 10% 
or less. 

To obtain the effective recombination coefficient for 
a certain line one has to weigh the upper level by the 
corresponding radiative ratefl As an example, we calcu- 
late the normalised effective recombination coefficient /eis? 
of the quinted 4d('^D) 3p(^P) 6157 A line at lOOOK. 
The radiative rate from 4d(^D) to 3p(^P) is given by 
^6157 = 7.62(-f6) s"^ and the radiative rate 4d(^D) to 
4p(^P) is given by j42651i = 6.44(-|-6) s~^ and therefore 

/6157 = /4d(SD)-i = 3.8(-2) (21) 

A(il57 + ^26511 

which agrees well with the value given by lJulienne et al.l 
(|1974 ) at 1160K. 



http://www.nist.gov/index.html 
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Table 1. Effective normalised recombination fraction of O I for triplet and quin ted states of O I n = 3. The 
rates are normalised on the total radiative recombination rate of lBadneOI l|2006l V States indicated by * show 
the case of optical thin ground state transitions. Values in brackets indicate a power of 10, e.g. a(b) = a ■ 10*. 
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